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ABSTRACT 


In this study, a self-consistent derivation of the conservation laws is given for flows of 
a fluid-solid mixture. A unified analytical framework for obtaining constitutive relations i 
nrovided This analysis uses a control volume/control surface approach that is widely use 
fn fluid mechanics. All terms in the governing equations and the 

written in terms of the mass-weighted averages except solid concentration. It is belie ed 
strongly that the mass- weighted average is the natural bridge b ?^en mici rom< echa nics d 

constitutive relations. The derived momentum equations contain terms that ditter from all 
existing models except that of Prospered and Jones (1984). However, their wumtou 
are nof needed here. P Special attention is given to the solid phase pressure. The physical 
basis of prevteusrty assumed form for this pressure (Givler 1987) becomes clear. A number 
of related phenomena are also discussed. These include the anti-diffusion and anisotropic 
normal stresses. The energy equations are also different from existing models. But detail 
discussion on the energy equations is left to future work. 


I. INTRODUCTION 

Modeling a flowing fluid-solid mixture starts from writing down a set of governing 
equations. These equations describe the conservation of mass, momentum and energy. In 
the early stage, the popular approach was to view the mixture as a single phase mate- 
rial. Consequently, the following type of equations were used for mass and momentum 
conservations (see Zuber 1964, Ishii 1975, 1977). 

f%l + V-(pm>u m = 0 (!) 

at 

4r(Pm)Um + V • (Pm)u m U m = V • T"* + (Pm) g ( 2 ) 

at 

In which, ( ) represents the ensemble average, p m is the mixture density, u m is the mixture 
velocity, T m is the mixture stress and g is the body force per unit mass. 
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Recently, the two-phase flow approach becomes more popular. Using this approach, 
the conservation equations are formulated for each individual phase separately. This ap- 
proach allows for direct modeling of a fluid-solid mixture when the two phases have dis- 
tinctly different dynamics. For mass and momentum, these conservation equations are 
generically written in the following way (see Ishii 1975, Drew 1976, 1983, McTigue et al. 

1986). 

d(p 3 c)/dt + V • (p,cu) = 0 , . 

d(pf(l-c))/dt + V'(p f (l- c )v) = 0 (4) 

d{p,cu)/dt + V • (p,cuu) = ( Pt cg) + (m) + V • (cT.) 

d{ Pf (l - c)\)/dt + V • ( Pf ( l - c )vv) = ( Pf ( l _ c ) g) - (m) + V • ((1 - c)Tf) (6) 

where p a and p f are densities of the solid and fluid phases, u and v are the velocities of 
the solid and fluid phases, m is the interaction force per unit volume of the mixture, c is 
the local solid concentration (equals to 1 at a solid point and 0 at a fluid point), T. and 
T f are the solid and fluid stress respectively. The above equations do not consider the 
phase changes at the interface. The energy conservation equations have not been studied 
as extensively as the mass and momentum conservation equations. 

Because the two phases are separated, available information on a single particle’s mo- 
tion and the particle-particle interactions in a fluid environment are incorporated directly. 

In a mixture model, these informations will first be utilized to obtain the drift flux term in 
T of Eq. (2). Preference of the two approaches apparently depends on the application. 
However, since the mixture model can be derived from adding together the two-phase 
equations, the two-phase approach is considered more fundamental. 

Many mathematical models have been derived based on Eqs. (3)-(6). There are two 
issues in modeling terms appeared in those conservation equations. First, what kind of 
averaging is used in defining the macro quantities, such as concentration, velocity, velocity 
divergence, strain-rate, etc. Second, how to obtain the required constitutive relations for 
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averaged terms like stress, pressure, phase interaction, etc. Neither of these two issues rs 
settled at the moment. 

We will refer to the mechanics that governs particle’s motion as the “micromechanics 
Such micromechanics includes particle-particle and particle-fluid interactions. In order to 
explicitly formulate constitutive relations, knowledge of micromechanics is essential. 

It is understood that the formulation of constitutive relation is extremely compli- 
cated, because the micromechanics of fluid-solid interaction is itself not well understood. 
Another reason for this difficulty is less obvious but more significant. That is, up to the 
moment, there has not been a set of governing equations in which all terms are interpreted 
with a consistent averaging method, without such governing equations, and a consistent 
bookkeeping to account for the micromechanics. It is impossible to correctly formulate the 
constitutive relations. The authors believe that this is the source of the recent argument 
about the “solid phase pressure” and related phenomena. A survey of the recent literature 
shows that different interpretations have been given to the solid phase pressure. It has 
been suggested to be equal to the (i) fluid phase pressure (Drew 1976), (ii) averaged fluid 
pressure around the surface of a particle (Givler 1987), or (iii) a more sophisticated version 
of (ii) with additional consideration of Brownian forces and bulk viscosity (McTigue et al. 
1986). All of which are intended to apply to an arbitrary flow of a fluid-solid mixture. 
Similarly, the phase interaction (m) in Eqs. (5) and (6) has also been modeled in many 
different ways. Essentially, in all the more recent works this term has been modeled as 

(m) = nh + pV(c) ( 7 ) 

where n is the averaged number of particles per unit volume of the mixture, h is the 
averaged hydrodynamic force per particle. The term p is a source of confusion again. It 
has been equated to the (i) averaged fluid pressure over the fluid-solid interface (Drew 
1983), (ii) fluid phase pressure (McTigue et al. 1986) or (iii) hydrostatic fluid pressure 
(Ahmadi 1987). Again, all the above are suggested for general flow condition. 
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Although the existing two-phase flow models appear to be inconsistent, development 

of the above models has provided a great deal of insights. These insights are essential to 
the work presented here. 

In the following, we will give a step by step derivation of the two-phase flow governing 
equations based on the mass-weighted average defined 


as 


\ C *) 


( 8 ) 


where c k equals to c for the solid phase and 1 - c for the fluid phase. This averaging 
method has been exploited in the theory of compressible turbulence (Van Driest 1951), 
and applied to granular flow (Ahmadi and Shahinpoor 1983). 

The mass-weighted average of any given quantity is the quantity averaged within that 
phase only. For instance, mass-weighted solid velocity at a point is the average velocity of 
all observed particles passing that point. This average is the easiest one to measure in most 
real flows. Moreover, through using this average, a direct bridge between micromechanics 
and the constitutive relations may be established. Applying this average to the solid phase 
stress, terms such as the solid phase pressure will have a clear meaning. Thus a unique 
definition for these quantities in terms of micromechanics is possible. In addition, a number 
of interesting phenomena that seem to defy the well accepted property of fluids are observed 
when a fluid-solid mixture is viewed as a composition of two separate continuums. 

In the present study, we concentrate on the derivation of governing equations, in- 
terpretation of the averaged terms and the development of constitutive relations from 
micromechanics. Derivation of the actual constitutive relations for various flow conditions 
is beyond our scope and present ability. A few exceptionally simple cases will however be 
studied. In order to complete the mathematical modeling, boundary conditions must be 
derived. Due to the existence of the discrete solid phase, derivation of boundary conditions 
is equally difficult as constitutive relations. This is also left to future development. 
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II. governing equations 


Consider an arbitrary control volume V as shown in Fig. 1. Its surface S is the 
control surface. For visual reason this control volume is drawn large compared with the 
solid particle's size. The derivation is not restricted to this size. The particles me not 
necessarily spherical or uniform either. However, in order to simplify the notations, we 
will discuss the case of uniform spherical particles. The essence of the analysis is captured in 
this simplified case. There has been previous work deriving conservation equations for this 
type of control volume for a Auid-solid mixture (Soo 1981). Nevertheless, such equations 
have not been derived in terms of the mass-weighted average nor the interpretation in 
terms of micromechanics given. 



Fig. 1 A control volume V with control surface 5 = 5, + 5/ and 
internal interface S{. 

In general, the flow situation is such that the condition in any given control volume 
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from different observations is a random phenomenon. Only an average behavior may be 
described by a deterministic mathematical model. To allow for cases where flow conditions 
may vary over a length scale comparable with the particle size, the ensemble average over 
many realizations of a control volume smaller than the “representative volume” commonly 
used in fluid mechanics may be necessary. 

For one realization, the rate of increase of solid mass in the control volume is 

^] v f>.cdv + j s e. m .t>ds = o (9) 

»here A is the unit outward normal. After ensemble averaging, the resulting integrands 
are smooth functions. Applying Green’s theorem to the second term and removing the 
integrals, the mass conservation of the solid phase in this control volume is obtained as 

d(p s ) 

~dt +*■{,•>{»} «0 ( 10 ) 


In the above, (p.c) is replaced by (p-) and (p.cu) is replaced by (p-){u} where { } is the 
mass-weighted average defined in Eq. (8). Similarly, for the fluid phase, 


8{pf) 

dt 


+ V - (p / ){v) =0. 


( 11 ) 


The momentum conservation equations require representation of the forcing terms 

which include the surface and body forces. This conservation for the solid and fluid phases 
in one realization are, respectively, 

§{ v P.™*V + J s p.cuu-nJS = J s cT..fidS + f v p.cgdV + J m JV (12) 

dt X P/(1 - c)v dV + J s M 1 ~ c ) vv ndS = j(\ - c)Tf ■ A dS 

+ f Ml- c)gdV- [ mdV (13) 


40 



Again the ensemble average is first applied to smooth the above integrands. The 
Green’s theorem may then be used to change the surface integral to volume integral. After 
removing the integral sign, the above equations become 

— (p 3 cu) + V • (p,cuu) = V • (cT.) + {p»cg) + <m) (14) 

dt 

1 - c)v) + V • (p f ( 1 - c)vv) = V • ((1 - c)Tf) + (Pf( 1 - C) g) - (m) (15) 

dt 

Substituting u by {u} + u" and v by {v} + v" and making use of the mass conservation 
equations, we obtain the following equations with the Reynolds stresses for both phases. 
In these equations mass- weighted velocities and stresses appear. 

(p 9 ){d{n}/dt + {u} • V{u}) = V • (c){T 8 ) - v • </){u"u"} + (p 9 ) g + (m) (16) 

(p f )(d{v}/dt + {v} • V{v}) = V • <1 - c){T f } - V • (p / ){v"v"} + (pf) g - (m) (17) 

The Reynolds stress in the solid phase is also called the kinetic stress in the granular flow 
terminology. 

In a realization over a period of time, particles cut by the control surface may in- 
teract with neighboring particles through collisions. The rate of momentum transfer to 
the interior of the control volume resulted from these collisions is part of the surface force 
T*. Moreover, the hydrodynamic forces acting on particles cut by the control surface also 
produce surface forces that contribute to T s . We denote these two stresses as T c and TP 
respectively. 

Concept for modeling the mass- weighted average of T c has been described in the 
granular flow literature (e.g. Bagnold 1954, Jenkins and Savage 1983, an excellent survey 
to appear by Campbell 1990). Although most of the work deals with negligible fluid 
effect, the route to extend to fluid-solid mixture is, though complicated, quite clear (e.g. 
Ackermann and Shen 1982, Shen et al. 1988). On the other hand, the explanation of 
the hydrodynamic stress on the solid phase is not readily available. Intuitively, one would 
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quite comfortably accept that the fluid only acts on each particle in the control volume 
through the drag, added mass, etc. Hence, it should contribute to the body force only. 
Consequently the effect of fluid-solid interaction is the hydrodynamic force per particle 
multiply the number of particles per unit mixture volume. This concept is proved doubtful 
as demonstrated by the existing various models discussed in Eq. (7). In the following we 
will rigorously formulate the fluid effect in a two-phase flow. 

Consider a surface particle P in Fig. 2. Part of this particle, P\ is inside the control 
volume, part of it, P° is outside of it. The hydrodynamic force acting on P produces a 
pair of internal forces, ± t, on the intersection of the particle and the control surface S a . 
The total hydrodynamic force acting on this particle, h, can similarly be decomposed into 
that on the outside of P, d°, and that on the inside of P, d'. 



Fig. 2 Decomposing the hydrodynamic force on a surface particle P. 

The total hydrodynamic force i6 h = d 5 + d°. 

In budgeting the total force acting on the solid portion of V, the interface force t 
acting on 5, most naturally belongs to the solid phase stress. In fact it may equally well 
be classified as part of the body force since after all it acts on P* which is inside of V. 
As long as the budgeting of all forces is done in a consistent way, the resulting equation 
of motion should not depend on the detail of the bookkeeping. Most of the existing 
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two-phase models appear to fail in using a consistent bookkeeping. That is, force acting 
on the control surface and that acting on the internal interface are not always carefully 
distinguished. We choose to call t the surface force for V and d' the body force for V due 
to fluid-solid interaction. 

The internal force t is the difference between d° and the total surface force acting on 
P°. The first may be determined if the hydrodynamic force distribution on a particle’s 
surface and the particle’s location relative to the control surface are both given. 

The second is determined using Newton’s second law. That is, the total surface force 
acting on the partial particle P° equals to its inertia subtract the body force acting on 
it. Ensemble average of the inertia and the body force on P° may also be obtained if 
the particle’s relative position with respect to the control surface is given. For a control 
volume reasonably away from the boundary of the flow field, the position of particles on 
the control surface S s may be assumed to uniformly distribute inside or outside the control 
volume. Using these arguments it has been shown in Hwang and Shen (1989a) that the 
ensemble average of the sum of t in a unit area produces the following stress, 

{TV] = vXh (18) 

where A 0 and V 0 are the surface area and volume of particle P respectively, E is the local 
fluid stress on the particle’s surface, and r is the position vector. The very same equation 
has been obtained by Batchelor (1970) with a quite different derivation. This term has 
been named the “particle-presence stress” in Hwang and Shen (1989a), and it contains the 
“interaction stress” in McTigue et al. (1986). The second term inside the parenthesis in 
Eq. (18) is related to the particle’s rotational inertia (Hwang and Shen 1989a). 

Next we discuss the interaction term m in the solid momentum equation. The most 
natural way consistent with the above bookkeeping is to define it as the total force acting 
on the interface of the two phases inside the control volume V. This interface is denoted 
by Si in Fig. 1, which consists of the surfaces of whole particles if they are entirely in V, 
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otherwise only the part of particle’s surface that is inside V belongs to 5,. Total force 
acting on 5, is the sum of all h, the hydrodynamic force on a whole particle, and all d‘, 
the portion of the hydrodynamic force on P inside of V. As shown in Hwang and Shen 
(1988), the ensemble average of the total interaction force m is 

(m) = v; {h} + (c)v ' {Tf} " v • C < TP } (19) 

In the above, {T f } = {T f } - is used to indicate the total fluid stress. Substi- 

tuting Eq. (19) and T. = T c -f- T p into Eq. (16) , the solid momentum equation becomes 

<p')(d{u}/c* + {u}.V{u}) = V.((c){T c })-V.(/){u"u"} + (/)g+^{h} + (c)V.{T f } 

( 20 ) 

This equation is identical to that derived in Prosperetti and Jones (1984). In that work, 
a couple of assumptions were made to arrive at Eq. (20). It is shown here that those 
assumptions are not necessary. Similarly, the fluid momentum equation may be derived as 

(/K9(v}/Sf+{v).V{v}) = (l-2(c))V.{T f }-{T f }V(c)+(/)g-M {h}+v .( (c){TP)) 

V Q 

( 21 ) 

As it shows in above equations, the effect of particle-presence stress {T p } appears only in 
the fluid momentum equation. 


Derivation of the energy equations follows the same spirit as the above. Hwang and 
Shen (1989c) gave the following equations in indicial form for the turbulence energy in the 
solid and fluid phases respectively, 

(p,) (£ {i hr~ } + {u} ■ v{iL ir" } ) = (c){TC} : v{u > + < c H TP > : v m 

- <P*){u"u"} : V{u} + V • ((c){T c • u"}) + V - ((c){T p • u"}) 

, u W U W 


- V • •»“}) + (». u") - < £ ){y } 


(22) 


(}L T L} + (v) ■ V <- !L T L}) = ~ (1 - c )(p / ) v - {v} - (/){v"v") : V{v) 

II ti tf 


-v((/){(^ + ^K} ) + ^« i-cH 
-M 1 - c){Vv" : Vv") - (m • v") 


Y ,; . 


}) 


( 23 ) 
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where 7 * is the sink of solid turbulence energy from particle-particle interaction (e.g. 
collisions! dissipation), p f is the fluid pressure, p f " = p f - { p f }, and p is the dynamic 
viscosity of the fluid. 

The interaction force m does unequal amount of work to the turbulence energy in the 
fluid and solid phases. This difference, as given in Hwang and Shen (1989c), comes from 
the slip velocity between the two phases. 

The combination of mass, momentum and energy conservation equations given by the 
six equations: Eqs. (10), (11), (20), (21), (22) and (23) governs the six unknowns c, u, v, 
{u" • u"/2}, {v" • v"/ 2} and { p f}. To form a closed system, a large number of constitutive 
relations must be obtained. Most of these constitutive relations require knowledge beyond 
the current understanding of fluid mechanics around a particle. One example is the term 
{m • u"} appeared in Eq. (22). As shown in Eq. (19), m contains the hydrodynamic force 
h acting on a particle. If one considers the drag force part of h, which is in general in 
terms of u" — v", the term {m • u"} will produce {v" • u"}. Namely, the correlation of 
the fluid and solid turbulence must be formulated as part of the constitutive relations. A 
number of analyses have indicated that this correlation depends greatly on the particle’s 
size and density (Chao 1964, Xie 1987, Abou-Arab and Roco 1988). Quantitative study 
of such correlation is far from complete. 

Despite the complexity and the lack of necessary knowledge in fluid mechanics, a 
number of results may be obtained for the constitutive relations appeared in the governing 
equations. In the next section, we will derive a term associated with the fluid-solid in- 
teractions. Namely, the particle-presence stress T p including a detailed discussion on the 
solid phase pressure. A number of interesting observations will be given at the end of the 
derivation. 
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III. CONSTITUTIVE MODELING 


To simplify the notation, we remove all the average symbols ({ ) and { }) and adopt 
indicia] notation from this section on. 

In this section, we will discuss the particle-presence stress Tf-. Detail of this term is 
given in Eq. (18). The term in Eq. (18) related to the particle’s rotational inertia cancels 
with the rotational contribution in the particle’s Reynolds stress in Eq. (16) (Babic 1989). 
Hence the net particle- presence stress is only the first term on the right side of Eq. (18). 
This term may be quantified if the distribution of the hydrodynamic stress E is known. 
Before quantifying this stress for special cases, we will provide an interpretation of the 
solid phase pressure first. 

As defined in Continuum Mechanics, pressure means the negative average of the nor- 
mal stresses. Therefore, in addition to the contributions from the collisional and Reynolds 
stress, the particle-presence stress will also produce the following solid phase pressure. 


pP = 


R 3 


3V C 


$ ?2n + ir 

- I / Y<iknkni sin <f>d<f>dd 

o Jo Jo 


l r** r n 2 

~ - 2 ^ e,l ^ ik +%l J ‘ e ik)nkni sin <f>d<t>d8 

= - f r 

4tt Jo Jo 


p J sin (f> d<f> dO 


(24) 


where 4> and 0 are the polar and azimuthal angles respectively. In the above, a Newtonian 
fluid is assumed such that 


Si* = ~P f 6ik - -fieuSi k + 2pe ik (25) 

where e** is the component of local strain-rate. Moreover, for Eq. (24) to be valid, the 
fluid flow near the particle must be incompressible and the strain-rate must possess certain 
symmetry property. Thus the viscous contribution vanishes from Eq. (24). Under these 
special conditions, it is seen that the solid phase pressure is the numerical average of the 
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fluid pressure around the particles. This very result has been postulated by Givler (1987) 
and derived by Prosperetti and Jones (1984) using a couple of assumptions. 

We now derive the particle-presence stress in two extreme cases for a very dilute 
mixture. The two extreme cases are: Stokes regime and inviscid regime. In both cases we 
assume the fluid is Newtonian and the flow is incompressible. 

In case of a shear flow in the Stokes regime, the local strain-rate around a spherical 
particle in an infinite flow field is (Batchelor 1967) 

5 R 3 5 R 5 \ 


, s R 5 \ , r, (rir k 8ji + rjr k Su 2r*r, ^ ( 5R 

=«<r0 = M 1 - 7 ?) + M 3~ e, 0 ( ■ 2 ^ + 


r 5 J 


+ 


( r ' r i 1 c \ 

/25 R 3 

35R 5 \ 

(tT - 3 6 -’) 

V 2r 3 

2r 5 ) 


(26) 


where rj is the ith component of a position vector r, r =| r |, Eij is the component the 


undisturbed strain-rate and R is the particle’s radius. 

The fluid pressure around a spherical particle in a uniform incoming flow Uoo is 
(Chester and Breach 1969) 

p = p> + ^ [-5( 1 + ) cos 4, + R. cos 2 * + 0(R\ log /?.)], (27) 

r z o uz 

where p f is the undisturbed fluid pressure; Uoo is equivalent to the relative velocity of the 
particle to the fluid flow, or Uoo = v r where v r =| v - u |; R t is the particle Reynolds 
number defined as 2 pfRv r /p] and <f> is the polar angle of a point on the surface of the 
particle measured from the direction of Uoo- 

Substituting Eqs. (25) and (26) into Eq. (18), the viscous contribution of the particle- 
presence stress is obtained as 

Tff = |(2p Eij) (28) 

As discussed in Batchelor (1970), this stress together with the viscous stress in the fluid 
phase reproduces Einstein’s formula for the effective viscosity (1956). The solid phase 
pressure is shown below after substituting Eq. (27) into Eq. (24), 

I?=P f + ^eivl ( 29 ) 
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In the above, we used •p p to indicate the fact that only particle-presence stress T p is 
included in the total solid stress. Neither collisions! nor Reynolds pressure is considered. 

The above is not the whole story about the particle-presence stress. If one by-passes 
the solid phase pressure and investigate the normal stresses directly, the following is found. 


m = v - ^ 

'TPP _ rrPP 
- l 22 ~ 


33 

/ 27 2 

= ' 160 ^ fVr 


(30) 


The total particle-presence stress is the sum of Eqs. (28) and (30). It is apparent that 
Eq. (29) can be reproduced by taking the negative average of Tff, T 2 P 2 P and T pp . 

Eq. (30) shows that on top of the viscous effect the isotropic fluid pressure induces 
anisotropic normal stress in the solid phase stress. This anisotropy of the normal stress is 
the product of the distinctly different dynamics of the two phases involved. If both phases 
move with exactly same velocity, this phenomenon will disappear. 

In case of an inviscid flow, only the fluid pressure contributes to the particle-presence 
stress. This pressure around the spherical particle is (Lamb 1932), 


1 19 

Poo + jP/t>r(-g + g COS 24,). 


(31) 


Substituting Eq. (31) into Eq. (25) and neglecting the viscous part, Eq. (18) becomes 

T’n = V - 

rrP rpp 

J 22 “ ** 33 

f 2 , 

= -p + ( 32 ) 


and the solid phase pressure from the hydrodynamic effect alone is the negative average 
of the above, 

P P = P f - \pfv 2 r (33) 
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The anisotropy of the normal stresses is again observed. Theoretical speaking, without 
considering Reynolds stress, one may observe the shear force in a two-phase mixture even 
though the fluid is inviscid. 

Comparing the coefficient of the solid phase pressure in the above two extreme cases, 
one would expect that as the particle Reynolds number increases from nearly zero to 
nearly infinity, the coefficient in front of the solid phase pressure should vary gradually 
from ^ to Explicit determination of this coefficient depends on the knowledge of the 
hydrodynamic stress distribution around the particle. Unfortunately this information is 
not available for the entire range of the particle Reynolds number. However, for particle 
Reynolds number greater than 3000, experimental data shows that the viscous contribution 
is negligible. In this case, sufficient information exists to empirically determine the pressure 
distribution, and accordingly the solid phase pressure (Hwang and Shen 1989b). 

IV. SOME INTRIGUING POINTS 


In deriving the governing equations, a few points have struck the authors as being 
quite non-trivial. Some of those may have significant implications that is unclear at the 
moment. 

First, the governing equations are for the mass-weighted quantities. All constitutive 
relations must eventually be written in these quantities to produce a closed system. The 
kinematic quantities appeared in these governing equations are the mass- weighted veloci- 
ties. On the other hand, in constitutive relations, we look for mathematical description of 
the stresses. The average stresses, at least for the viscous component of the fluid stress, 
is however a function of the average local fluid strain-rate instead of the gradient of the 
mass-weighted velocity. Namely, this stress is in terms of the mass- weighted strain-rate in 
stead of the gradient of the mass-weighted velocity. For a fluid-solid mixture with rigid 
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particles, this relation has been derived as (Hwang and Shen 1989c) 

,dvj 9{uj} 9(1 -c) {v.} 19, , / c \ 

dxi — dxi + dxi (1 - c ) + (1 - c ) 9i, “ (TrjjW'j} ( 34 ) 

where {fljj} is the spin tensor (or particle rotation) of the solid phase. In addition to 
the above two possible definitions for the phase strain-rate, a third strain-rate is also 
important. This is called the bulk (or mixture) strain-rate defined by 


E H = (1 - c){e f ij} + <c){e a 0 } 

= (1 - c){e { ij} (35) 

where and e*,j are the local strain-rate for the fluid and solid phase respectively. In 
case of rigid particles, e^=0. In representing the effective viscous stress in a fluid-solid 
mixture, the above bulk strain-rate must be used. 

In a pure fluid flow, such distinction shown in Eq. (34) or Eq. (35) does not exist. In 
a two- phase flow, this point becomes essential in many places. In fact, if one interprets 
the bulk viscous stress using d{vj}/dx x instead of {dvj/dxi}, the coefficient of effective 
viscosity in Eq. (28), i.e. |, will become § or | depends on different interpretations of the 
mixture strain-rate. 

Second, in Givler (1987), a phenomenon called anti-diffusion in the solid phase was 
mentioned. This term is believed to come from an inconsistent bookkeeping system of the 
surface and interaction forces in the control volume V. The equation given there for the 
solid phase was 


+ “ii)) = + yK - eg - (pP - pf)± 


(36) 


When the potential flow is considered, the 9c/9x, term has a positive coefficient. 
This presents a force that moves the solid phase from a low concentration towards high 
concentration, which is a puzzling phenomenon. However, in the present derivation, it 
is shown that in the solid momentum equation such dc/dx x term does not exist. The 
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associated diffusive or anti-diffusive forces are thus absent. But, investigating the fluid 
momentum equation, we do observe such dc/dxi term. Through the mass conservation, 
any diffusion in the fluid phase results an effective anti-diffusion in the solid phase and vice 
versa. At this point, we have observed a totally opposite trend between Givler (1987) and 
the present model. Namely, Givler’s model would predict a diffusion phenomenon for solid 
particles in the Stokes regime and anti-diffusion phenomenon in inviscid flow regime. The 
present model however predicts exactly the opposite. Further verification with detailed 
experimental data is desirable. 

On top of the above, the anisotropy of the solid phase normal stress produced by fluid 
pressure is a new observation. This may produce interesting thermodynamic properties 
that axe peculiar to a two-phase flow only. 

It is natural to ask whether such a great care in deriving governing equations shown 
here is of importance. In order to see this, two models have been applied to a vertical pipe 
flow of a fluid-solid mixture (Hwang and Shen 1988). The two models are identical except 
the phase interaction term m,. For case (A), rrii is modeled as shown in Eq. (7) with p = p p 
and case (B) is the present model given by Eq. (19). The resulting non-dimensional slip 
velocity u * — v* verses the non-dimensional fluid pressure gradient (^f )* is reproduced in 
Fig. 3 for three different density ratios p* = ^ . Qualitatively different results are obtained 
from these two models. From the behavior of the neutrally buoyant particles, it is believed 
that the present model is more reasonable. 

V. CONCLUSION 

For a two-phase flow, the interpretation of terms in the governing equations is the first 
step towards deriving constitutive relations using micromechanics. Therefore it has been 
surprising to us that the work presented here has not been available in the vast amount of 
two-phase flow literature. On the other hand, this may be explained by the fact that only 
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recently researchers have started to derive the constitutive relations from micromechanics. 
Hence the need for such derivation of governing equations is also very recent. 



Fig. 3 Effects of p, m and (dp/dx) m on slip velocity from 
model A and model B. 

It is shown in this work that the derivation of governing equations for flows of a fluid- 
solid mixture is not as straightforward as any single phase continuum. This derivation 
is done with a careful and almost tedious method to account for all transfer quantities 
between each phase. The resulting equations in terms of mass-weighted average quantities 
differ from the existing models. 

These governing equations provide a foundation to incorporate micromechanics in 
deriving the constitutive relations. A number of previously well accepted facts about fluid- 
solid flows emerge naturally from this analysis. In addition, interesting phenomena such 
as: non-equality of average strain-rate and gradient of the average velocity, anti-diffusion, 
and anisotropic normal stress have been observed in these equations. Interpretation and 
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implications of these phenomena is of great interest. 
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